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Abstract. - The in- and out-of- plane lower critical fields and magnetic penetration depths of 

(-H ■ LiFeAs were examined. The anisotropy ratio Th^^ (0) was smaller than the expected theoretical 

Q ' value, and increased slightly with increasing temperature from 0.6Tc to T^. The small degree of 

^ , anisotropy was numerically confirmed by considering electron correlation effect. The temperature 

dependence of the penetration depths followed a power law (~T") below O.STc, with n>3.5 for 

both \ab and Ac. Based on theoretical studies of iron-based superconductors, these results suggest 

that the superconductivity of LiFeAs can be represented by an extended s±-wave due to the weak 

v^ ' impurity scattering effect. And the magnitudes of the two gaps were also evaluated by fitting the 

f^ ' superfiuid density for both the in- and out-of-plane to the two-gap model. The estimated values for 

0\ . the two gaps are consistent with the results from angle resolved photoemission spectroscopy and 

^^ ' specific heat experiments. 

o 
o 

Introduction. — Since the discovery of iron-based superconductors, considerable ef- 
l^% ' fort has been devoted to identifying a higher transition temperature (Tc) and examining 

^ . the mechanism for the superconducting (SC) behavior. Determining the symmetry of the 

superconducting order-parameter (OP) is an important first step towards unveiling the SC 
mechanism in iron-based materials. The magnetic penetration depth X{T) is one of the most 
useful parameters for probing the symmetry of the OP. The SC gap structure and symme- 
try for 1111- and 122-type superconductors have been studied by estimating the magnetic 
penetration depth using tunneling diode resonators (TDR) [lH3] , microwave methods [IHB] , 
and scanning SQUID susceptometry [71[S]. Previous studies on the magnetic penetration 
depth in SmFeAsOi.8Fo.2 0, PrFeAsOi_y [S], and Bai_2;K2;Fe2As2 [5] reported an expo- 
nential temperature dependence, supporting the existence of a s-wave-type SC gap. On the 
other hand, the magnetic penetration depth for LaFePO [7] exhibited a near linear temper- 
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ature dependence, supporting a nodal gap. Recent reports on Fe-based superconductors, 
however, revealed a power law behavior with an exponent of 2<n<2.5 at low temperatures 
[Hin]. Nevertheless, symmetric properties of SC gaps, such as extended s-wave [5], s±-wave 
[T2HI3, nodal d-wave [TIFH]. and point nodal gap [3], have been suggested. Recent theoret- 
ical advances |151ll6| that consider an impurity scattering effect in the Eliashberg equation 
have provided convincing evidence that the s-t-wave SC state is the most promising can- 
didate for the true pairing state of the Fe-pnictide superconductors. For anisotropy of Fe 
based superconductors, experimental measurements have revealed very small anisotropy ra- 
tios compared to those of cuprate, although both superconductors have the same layered 
structure [T71[TH]. Furthermore, some theoretical studies have suggested the importance of 
electron correlations [19] and their effects on the unconventional superconductivity [20l[2T] 
of Fe pnictide. 

In this letter, we show that both magnetic penetration depths, Xab for Ii\\c and Ac for 
H\\ah, follow a power law (~T") with n>3.5 below CST^. These results suggest that the 
superconductivity of LiFeAs can be represented by an extended s±-wave model due to a weak 
impurity scattering effect. In addition, the anisotropy ratio of lower critical field is found to 
be smaller than the expected theoretical value. According to first principles calculations for 
LiFeAs compounds, the anisotropy ratios of rA(0) and rp(0), estimated using the magnetic 
penetration depth and electrical resistivity are approximately 3 and 9, respectively [22] . We 
also show that the anisotropy ratio can be examined by considering the electron correlation 
effect on the band structure using a local density approximation (LDA) and dynamical 
mean field theory (DMFT) approach. Finally, we report that LiFeAs possess two SC gaps 
determined through fittings to superfluid density for the in- and out-of-plane directions. 

Experiment. — The single crystal growth of LiFeAs with T^'^'^"—!! .b K is detailed 
Ref. [53] . The single crystal used in this study was obtained from another batch of the same 
ingot referred to in Ref. [23]. The magnetic field dependence of the magnetization M[H) 
was scanned from -70 to 70 kOe every 0.5 K from 2 to 7 K, and every 1 K from 7 K to Tc, 
for both H\\ab and H\\c using a vibrating sample magnetometer (VSM SQUID, Quantum 
Design) . 

Results and Discussions. — The inset in Fig. 1(a) shows the field dependence of the 
local magnetization measured at 2 K after zero field cooling. The shape of the magnetiza- 
tion loop is almost symmetric about the horizontal axis, which suggests that the hysteresis 
arises mainly from bulk flux pinning rather than from the surface barrier. Fig. 1(a) shows 
magnetization curves M{H) for H\\ab in the low field range below Tc- The M{H) data for 
H\\c shows similar behavior(not shown here). The value of H*i was obtained by evaluating 
the deviation of magnetic field from the Meissner line in the initial slope of the M{H) curve. 
As shown in Fig. 1(a), the magnetization very gradually deviates from the Meissner line, 
making the accurate location of H*i difficult to locate. To avoid any ambiguity in locating 
H*i from these broadened M-H curves, the following expression is employed 

Bcg(il)- (2A)i/2(iy + 4^Af)i/2 {H,i<H<H*). (1) 

Here, Beq{H) is the equilibrium flux density, yl is a constant, H* is the field defined as 
Beq{H*)=A, and H+4:TtM{= AttAM) is the amount corresponding to the deviation of the 
observed magnetization from the Meissner line. Eq.(l) indicates that, in principle, the plot 
of (AM)^/^ versus H should give the equilibrium B-H curve. Therefore, the lower critical 
field is given by the threshold field of this plot [53]. On the other hand, Eq.(l) assumes that 
the current density is field-independent. A previous study reported that the critical density 
of LiFeAs is independent of the magnetic field in low field regions below 700 Oe [23] , and this 
provides a reliable way of determining the value of H*^. Fig. 1(b) shows plots of AM versus 
H and (AM)^'^ versus H. The threshold of the (AM)^'^ plot may be slightly smeared due 
to the demagnetization field at the sharp corners of the specimen. Therefore, H*^ can be 
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Fig. 1: (Color online) (a) Initial region of tiie magnetization curve M{H) of single-crystal LiFeAs 
for H\\ab at various temperatures. Inset: M-H hysteresis loop at 2 K (b) AA/ as a function of H 
and (AA/)^' ^ as a function of -ff at 2 K, where AAf is the difference between the Meissner line and 
observed magnetization. 
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Fig. 2: (Color online). Extracted Hci of LiFeAs as a function of T for H\\ab and H\\c. The bars 
show the uncertainty of estimated by the deviating point of the M-H curve from the Meissner line in 
Fig. 1(a). The inset shows the temperature dependence of the anisotropy of Hci : ^H^i~Hrl/ H^'^ . 

determined by extrapolating the linear part of the (AM)^/^ versus H curve to intersect the 
horizontal axis i25j . The H^^ value is different from the actual lower critical field Hd due to 
geometric effects. The well-adapted relationship between H*-^ and Hd proposed by Brandt 
pS] was adopted in this case: Hci—H*i/i&vii[yjQ.'i(!)b/a), where a and h are the width and 
thickness of the slab-like superconductor, respectively. Here, a=2.8 mm and 6=1.3 mm for 



\\c_ 



if||c, whereas a=1.3 mm and &=2.8 mm for H\\ab. These values then yield H^^—2.b%H^ 



A\c* 



-\\ab_ 



r||ab: 



and Wd =l-41^ci j respectively. The lower critical fields parallel and perpendicular to 
the plane obtained for each temperature are plotted in Fig. 2 as a function of temperature. 
Hd (0) was evaluated by extrapolating Hd (T) in the low temperature region. As shown 
in Fig. 2, the flat behavior below 0.2Tc and the existence of an inflection point near 0.5Tc 
suggest that LiFeAs is a superconductor that has a full gap without nodes and a two-gap 



rll' 



Wab 



structure. The inset in Fig. 2 shows that the anisotropy of Hd , ^H^i =Hd I Hd j ^^^ ^ 
small value (between 1.2 and 2), and that Th^i increases slightly with increasing temperature 
above 0.6Tc. The temperature dependence of Th^i was similar to that observed for MgB2 
P7] and PrFeAsOi_y @1[S], indicating multiband superconductivity. Such behavior may be 
due to the existence of Fermi surfaces (FSs) of different sizes and anisotropics [35] . 
The following London formulae were used to estimate the penetration depth ; 



-\W_ 



A\ab 



-^ci— (*^o/47rA^j)lnKc, and Wd =('I'o/47rAa6Ac) luKaf,, where $o is the flux quantum and 
the Ginzburg-Landau parameter Kc=Xab/^ab, K,ab=\/ ^ab^c/ ^ab£,c- Figure 3 shows the tem- 
perature dependence of the magnetic penetration depths, Xab for H\\c and Ac for H\\ab, of 
LiFeAs. The depths at T=0 were obtained by extrapolating A(T) in the low temperature 
region, and were calculated to be Aa&(0)«198.4 nm and Ac(0)«250 nm. As shown in the 
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Fig. 3: (Color online) Temperature dependence of the magnetic penetration depths Xat and Ac 
of LiFeAs. Inset: Uemura plot of the superconducting transition temperature T^ versus \^^ for 
LiFeAs. 



inset of Fig. 3, although these are small errors associated with these values, the relationship 
between Tc and A~jf (0)=25.4 fira^^ is well represented by the electron-doped Uemura plot 
(dotted red line), and is also consistent with the results for the transverse-field muon-spin 
rotation (TF-/j,SR) of two polycrystalline LiFeAs compounds [33]. The anisotropy of the 
penetration depth rA=Ac/Aaf, exhibited small values of 1.26 at K and 1.8 near Tc;similar 
to the behavior seen for Th^i ■ 

To confirm the small anisotropy of the LiFeAs compound, the electron correlation effect 
on the band structure was examined using a LDA+DMFT approach. The LDA calcula- 
tions were performed using the full-potential linearized augmented plane-wave method in 
the Wien2k code [30] , and DMFT was implemented based on the LDA Hamiltonian method. 
The impurity problem within the DMFT self-consistency equation was solved using the con- 
tinuous time quantum Monte Carlo calculation. Parameters that give excellent agreement 
with the ARPES experiments [21] were used: C/=5.0 eV for the Coulomb interaction, and 



Table 1: Calculated renormahzation factors for each Fe 3d orbital obtained using the LAD+DMFT 
approach. 
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J=0.54 eV for the Hund's rule coupling constant. In addition, the specific heat coefficient 
7 was calculated to be 29.1 mJ/mol-K^. This is in good agreement with the experimental 
value of 35 mJ/mol-K^, obtained from C/T=j+/3T'^ in the measurement of the low temper- 
ature specific heat for the our sample (data not shown). Owing to the electron correlation 
effect, the ARPES of LiFeAs showed strong band renormalization: a factor of three larger 
compared to that of the LDA bands. Because each Fe 3d orbital has a different hybridiza- 
tion strength, each orbital also has a unique renormalization factor m* /m (where m* is the 
effective mass). Table 1 lists the calculated renormalization factors of each Fe 'id orbital. 
The smallest renormalization factors are those of the dz^ and d^^-y^ orbitals, whereas the 
dxy orbital has the largest renormalization factor. Owing to the difference in the renormal- 
ization factors, the band renormalization can be varied according to the orbital character 
in momentum space. Because the electrical anisotropy is strongly dependent on the band 
dispersion anisotropy, different band renormalizations can result in a large modification of 
the electrical anisotropy compared to that obtained from LDA calculations. Indeed, the 
LDA-t-DMFT calculation showed that pd Pab= 2.19, which is significantly smaller than the 
result of Pel Pab=^ obtained by LDA [35]. The result from the LDA-I-DMFT approach for 
LiFeAs corresponds with the electrical anisotropy ratio estimated from the experimental 
results: rp'-^3.3 [23] . 

Fig. 4 shows values of AA(T) estimated from the superfiuid densities p*(r)=iJci(r)/iJci(0)| 
=A^(0)/A2(T) of the in- and out-of-plane directions; note that AA(T)=A(T)-A(0). In the 
low-temperature region below 0.3Tc, the power law AA(r)~(r/rc)"' (solid blue line) pro- 
vides a suitable fit for both /^\ab{T) and AAc(T), with n=3.8±0.1 for H\\c and n=3.6±0.1 
for H\\ah, respectively. In previous reports on doped Fe-based superconductors, the penetra- 
tion depth behavior could be modeled using a power law with an exponent value of 2<n<2.5 
in the low temperature region [2l[^ . Many authors of theoretical articles have argued, based 
on the power-law behavior of the penetration depth at low temperatures, that most pnic- 
tides are s±-wave superconductors due to the impurity scattering effect [15j[T6]. According 
to Ref. [16] , the low temperature power- law behavior of the undoped LiFeAs superconduc- 
tor, with an exponent of n>3.5, suggests that the superconductor could be represented by 
an extended s±-wave model due to the weak impurity scattering effect. These results also 
correspond closely to the recent results of NMR and NQR experiments for polycrystalline 
Lia;FeAs compounds [35]. 

Recently, the s-wave two-gap feature has been examined using angle resolved photoemis- 
sion spectroscopy (ARPES) [31], specific heat [33], and lower critical field [36| experiments. 
In this study, the following expressions were used to examine the two-gap feature of LiFeAs: 

Ps{T)^xp,{T, Ai(0)) + (1 - x)p,{T, A2(0)) (2) 

p,{T)=l-2!-^^{-'4ll)^^,M==dE (3) 

Here, x (0<a::<l) is a weighting factor that indicates the contribution of the large gap, 
and f{E)=\\+exp{E/kBT)]^^ is the Fermi function. The temperature dependence of the 
gap A(r) can be expressed as A(r)=A(0)tanh[1.82[1.018(Tc/r-l)]°-5i] [34]. 

As shown in Fig. 5(a) and (b), the results of the fitting with two gaps (solid black 
line) were more consistent with the experimental data than those with one gap (dash- 
dotted blue line). The following parameters were used in the model: Aj''(0)=2.9±0.2 meV, 
Af (0)=1.4±0.1 meV, a;=0.55 for H\\c, and Af (0)=2.9±0.2 meV, A^(0)=1.2±0.1 meV, 
a;=0.51 for H\\ab. These results were in good agreement with those obtained using ARPES 
[3T1[35], specific heat [33], lower critical field [36], microwave surface impedance measure- 
ments [37], and TDR experiments [38]. Sasmal et al [36] reported the lower critical field, 
anisotropy and two-gap features of LiFeAs using a method similar to that reported in Ref. 
[55] . In the present study, however, the temperature dependence of Th^i due to the multi- 
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Fig. 4: (Color online) Temperature dependence of AA(T)/AA(0) for the penetration depths Xab 
and \c. The solid line(blue) is AA ~ {T /T^)" . 
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band effect, SC gap symmetry and isotropic properties of the gaps are novel, and the small 
measured anisotropy ratio was also confirmed using the LDA+DMFT approach by consid- 
ering the electron correlation effect. On the other hand, these results are in contrast with 
the single gap analysis reported by Inosov et al [3D]. 

Summary. — In summary, this study has confirmed the small anisotropy ratio of 
the lower critical fields and penetration depth in LiFeAs by both magnetic measurements 
and numerical calculations using the LDA+DMFT approach and considering the electron 
correlation effect. The low temperature power-law behavior (n>3.5) of the penetration 
depth suggests that the LiFeAs superconductor can be represented using an extended s±- 
wave model due to the weak impurity scattering effect. In addition, the sizes of the two SC 
gaps in LiFeAs for the in- and out-of-plane were obtained. The results were consistent with 
those reported previously. 
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